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ABSTRACT 

The Spitzer Survey of Stellar Structure in Galaxies (S^G) is a volume, magnitude, and size-limited 
survey of 2352 nearby galaxies with deep imaging at 3.6 and 4.5 /xm. In this paper we describe 
our surface photometry pipeline and showcase the associated data products that we have released 
to the community. We also identify the physical mechanisms leading to different levels of central 
stellar mass concentration for galaxies with the same total stellar mass. Finally, we derive the local 
stellar mass-size relation at 3.6/xm for galaxies of different morphologies. Our radial profiles reach 
stellar mass surface densities below ^ 1 Mq pc“^. Given the negligible impact of dust and the almost 
constant mass-to-light ratio at these wavelengths, these profiles constitute an accurate inventory of the 
radial distribution of stellar mass in nearby galaxies. From these profiles we have also derived global 
properties such as asymptotic magnitudes (and the corresponding stellar masses), isophotal sizes and 
shapes, and concentration indices. These and other data products from our various pipelines (science- 
ready mosaics, object masks, 2D image decompositions, and stellar mass maps), can be publicly 
accessed at IRSA (http://irsa.ipac.caltech.edu/data/SPITZER/S4G/). 
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1. INTRODUCTION 

Understanding how galaxies acquired their baryons 
over cosmic time is a key open question in extragalac- 
tic astronomy. How and when did galaxies of different 
types assemble the bulk of their stellar mass? Nearby 
galaxies are of particular relevance in this context: they 
constitute the present-day product of billions of years of 
evolution, so the past assembly history of these galax¬ 
ies is encoded in the present-day spatial distribution of 
old stars within them. Therefore, an accurate census of 
the current stellar structure of nearby galaxies provides 
essential constraints on the physics of galaxy formation 
and evolution. To address this critical issue, in this pa¬ 
per we present deep mid-infrared radial profiles for the 
more than 2300 galaxies in the Spitzer Survey of Stellar 
Structure in Galaxies (S'^G, see lSheth et al.ll20Tnl for the 
full survey description). 

Large surveys of nearby galaxies traditionally have 
been carried out in the optical regime, for example the 
Third Reference Catalogu e of Bright Galaxies (RC3, 
Ide Vaucouleurs et al.lll99 1fl. t he Sloan Digital Sky Sur- 
vev ISDSS. lYork et al.l 1200011. or the Carnegie-Irvine 
Galaxy Survey (|Ho et all 1201111 . However, translating 
optical measurements into stellar masses is not a straight¬ 
forward task. On the one hand, the optical mass-to-light 
ratio (M*/L) is strongly d ependent on the star f ormation 
history of the galaxy le.g. lBell fc de ,Ton^l2001[l . On the 
other hand, internal extinction by dust obscures a sig- 
nificant fracti on of the optical output of galaxies (e.g., 
lCalzettill2001ll . Observations at near- and mid-IR bands 
can circumvent these issues. At these long wavelengths, 
M*/L is only a shallow function of the star formation 
history, and dust extinction plays a minor role, leading 
to milde r variations in AA /L compared to optical wave¬ 
lengths (|Meidt et al.ll20l4l . 

Several wide field IR surveys have been carried out 
over the yea rs, such as the Two Micron All Sky Sur¬ 
vey ( 2 MAS S. lSkrutskie et al.ll2006l see also lJarrett et al.l 
20011 1200311. the De ep Near Infrared Survey (DENIS, 
Epchtein et al.l[l9^ and more re cently the Wide-fiel d 
Infrared Survey Explorer (WISE, iWright et al.l 1201011 . 
These surveys provide robust number statistics, but at 
the expense of a shallow image depth, insufficient to 
map the stellar content in the faint outskirts of galax¬ 
ies. In the case of WISE, the spatial resolution at 
3.4 and 4.6 p,m is relatively coarse (~ 6"). Conversely, 
other surveys have obtained much deeper IR images 
for small samples of several tens or a few hundred ob¬ 
jects. Projects like these include the O hio State Uni¬ 
versit y Bright Galaxy Survey lOSUBGS. lEskridge et al.l 
200211. the Near-Iiifrared atlas of SO-Sa galaxies (NIRSOS, 
Laurikainen et all 120111) . the Spitzer Infrared Nearby 
Galaxies Survey fSINGS. iKennicutt et al.ll2003ll and the 
Spitz er Local Volume Legacy Survey fLVL. iDale et al.l 
1200911 . But a complete understanding of galaxy assembly 
involves many independent parameters, such as galaxy 
mass, morphology, environment, bar presence, gas and 
dark matter content, etc. With so many independent 
dimensions along which a sample should be sliced, any 
sample with a few hundreds objects will be broken down 
into bins too small for a reliable understanding of the 
impact of any given parameter. 

S'^G was specifically designed to answer the need for 


a deep, large and uniform IR survey of nearby galax¬ 
ies. The survey contains over 2352 galaxies within 
40Mpc, away from the galactic plane (|6| > 30°), with 
an extinction-corrected B-band Vega magnitude brighter 
than 15.5 and a B-band diameter larger than 1'. We 
observed these galaxies at 3. 6 and 4.5 ixm wit h the In¬ 
frared A rray Camera (IRAQ , iFazio et al.ll2004ll onboard 
Spitzer (I Werner et al.ll2003l . We followed the success¬ 
ful observing strategy of the SINGS and LVL programs, 
reaching azimuthally-averaged stellar mass surface den¬ 
sities < IMqpc"^; this is a physical regime where the 
baryonic mass budget is dominated by atomic gas. 

The S'^G images are processed through a suite of differ¬ 
ent pipelines designed to produce a wealth of enhanced 
data products. The first three pipelines are described 
in this paper. Pipeline 1 (PI) creates science-ready mo¬ 
saics by combining all the individual exposures of each 
galaxy. Pipeline 2 (P2) masks out foreground stars, 
background galaxies and artifacts. Pipeline 3 (P3) per¬ 
forms surface photometry on the images and derives inte¬ 
grated quantities such as asy mptotic magnitu des, isopho- 
tal si zes, etc. Pipeline 4 lP4.[Salo et al.]l2015[ l uses GAL- 
FIT (|Peng et al.l I2002L1201011 to decompose each galaxy 
into bulges, disks, b ars, etc. Finally, Pipeline 5 (P5, 
iQuereieta et al.l 1201411 combines the 3.6 and 4.5/rm im¬ 
ages to produce stellar mass maps via an Independent 
Gompo nent Analysis , follow ing the methodology devel¬ 
oped by iMeidt et al.l (|2012ll . 

The goal of this paper is twofold. First, we describe the 
technical details and inner workings of PI, P2 and P3, as 
well as the resulting data products. Then, we discuss two 
particular scientific applications of our data: the local 
stellar mass-size relation, and the physics behind the vast 
diversity of morphologies in galaxies with the same stellar 
mass. 

For each galaxy we have obtained radial profiles both 
with fixed and free ellipticity and position angle. From 
these profiles we have also derived global measure¬ 
ments such as asymptotic magnitudes and stellar masses, 
isophotal sizes and ellipticities, and concentration in¬ 
dices. This dataset constitutes a unique tool to address 
many important issues on stellar structure, including but 
not limited to: 

• the radial structure of dark matter in galaxies. Our 
profiles constrain the radial distribution of stellar 
mass, which is a necessary ingredient when mod¬ 
elling rotation cu rves to infer the radial distribution 
of dark matter (iBosmal I1 978I iRu bin et al.l 119781 
iSofue fc Rnbinll^Oll: Ide Blok et al.ll2008|l. 

• the scaling laws of disks (iCourteau et al.l 1200711 , 
in pa rticular the loca l stel l ar mass-size rela¬ 
tion (|Kauffmann et al.l 1200,31: iShen et Ml 1200,311 
and, when kinematic d ata are availab le , the 
Tully-Fisher relation dTully fc Fisher! IT977t 
Aarqnson et al.l 1197^ iVerheiienI 120011 : iSorce et al.l 
20121: iZaritskv et al.ll2014ll! 

• the inside-out assembly of disks, by comparing 


our profiles of old stars with UV and opt ical pro¬ 
files probing younger stellar populations (Ide Jond 


19961 iBell & de Jonell2000l IMacArthur et al.l 

2004; 

Tavlor et al. 20051 Munoz-Mateos et al.l 

2007, 

2nnHWa,ng et aLIDnilll. 
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the local bar fraction an d the sizes, st re ngths 
and shapes of bar s ( Eskridge et al.l 


Erwinl 20051 IMenendez-De 

mestre et a 

Marinova & JoEfed |2(I(I7|: 

Sheth et al. 


2000 ; 


2002 ; 


2007; 


2ooa 


• radial migration of old sta rs and the a ssem- 


blv of galactic outskirts (IRoskar et al.l 

2008 

Sanchez-BlazQuez et al.ll2009l: IMinchev et al.l 

2011 

Martfn-Navarro et al.l 120121), in particular the 


study of breaks in radial projected surface den- 
sity profiles and their lin ks to bar resonances 
(|Munoz-Mateos et alJl2013l l. 


• one-dimensiona l _ str uctural decomp o sition s 

(jBaggett et al.l Il998t IMacArthur et alJ 1200311 . 
from which one can derive physical parameters 
such as disk scale-lengths, bulge effective radii, 
bulge-to-disk ratios, etc. 


the intrinsic face-on circularity of d isks a n d their 
vertical th ic kness (iSandage et al.l 119701: iRvden 


2004, 120061 iComeron et al 


ge et al.l liiJYUt It 
1.1 120111 iZaritskv' 


et al 


MU- 


• quantified galaxy morphology via non- 
parametr ic estimators such as concen t ration 


indices (iBershadv et al.l 120001: IConselicd 12003 

Abraham et al. 

120031: Munoz-Mateos et al.l 12009a 

Holwerda et al. 

fiiT1.l2014H. 


• the internal structu re of elliptical galaxies 
(iKormendv et ak ||2009[) , in p articular their intrin- 
sic axial ratios (iBvdenI 119^ . their Sersic indices 
(CaQn_eE5]|1 993ll and their boxy/disky structure 
( Bender et al.lll988l IPeletier et ai.lll990l i. 

This paper is structured as follows. In Section [2] we 
summarize our sample selection criteria and observa¬ 
tions. Section [3] contains a brief outline of the image 
reduction performed by PI. A more detailed descrip¬ 
tion can be found in Appendix [A] The object masking 
procedure is described in Section 0] Then, in Section [S] 
we explain our ellipse fitting technique and the result¬ 
ing data products. Readers interested in the scientific 
applications of these products can proceed to Section [ 6 l 
where we discuss the stellar-mass size relation, as well 
as the variety of morphologies and radial concentration 
that we find at a fixed stellar mass. Then, Section [7] ex¬ 
plains how to access and download our data. Finally, in 
Section [ 8 ] we summarize our main conclusions. 


against gas-poor early-type galaxies. However, our re¬ 
cently approved Cycle 10 program (prog. ID 10043, PI: 
K. Sheth) will complete Spitzer’s legacy with archival 
and new observations of all ^ 700 early-type galaxies 
within the S'^G volume, following the same observing 
strategy and selection cuts. 

Out of the 2352 S^G galaxies, ^ 25% were already 
present in the Spitzer archive as a result of previous pro¬ 
grams carried out during the cryogenic phase. We ob¬ 
served the remaining ^ 75% during the post-cryogenic 
mission. In this regard, from now on we will refer to the 
galaxies in our sample as either “archival” or “warm”, 
respectively. 

We observed the warm galaxies with a total on-source 
exposure time of 240 s per pixel, and mapped them out 
to at least 1.5 x D 25 . Depending on the apparent size 
of each galaxy, we used either a single dithered map or 
a mosaic of several pointings. Each galaxy was observed 
in two visits separated by at least 30 days, in order to 
image the galaxy with two different orientations thanks 
to the rotation of the telescope. This allowed for a better 
correction of cosmetic effects, cosmic ray and asteroid 
removal, and subpixel sampling. 

Most of the archival galaxies had been also mapped 
with at least 240 s per pixel and out to 1.5 x D 25 . How¬ 
ever, six archival galaxies had exposure times between 
90 and 200 s per pixel, and 125 were mapped out to less 
than 1.5 X D 25 . They represent a small fraction of the to¬ 
tal sample, and only specific science goals are affected by 
this. We therefore decided not to reobserve these galax¬ 
ies, as the incremental gain of repeating these observa¬ 
tions would not hav e made up for the required additional 
observing time ('see lSheth et al.l[2010l for more details on 
the sample selection). 

3. PIPELINE 1: SCIENCE-READY IMAGES 

We refer the reader to Appendix for a more de¬ 
tailed description of PI. Briefly speaking, PI creates 
science-ready mosaics by combining the different expo¬ 
sures of each galaxy. The pipeline first matches the 
background level of the individual exposures, using the 
overlapping regions among them. It then combines 
all frames following sta ndard dither/drizzle procedures 
(iFruchter fc Hookl2002fl . The final science-ready mosaics 
are delivered in units of MJysr“^, with a pixel scale of 
0.75". The FWHM of the PSF at 3.6 and 4.5 /rm is 1.7" 
and 1.6", respectively. For the farthest galaxies in our 
sample at ^ 40 Mpc, this translates into a physical size 
of ^ 300 pc. 

Besides the scientific images themselves, PI also pro¬ 
duces weight-maps showing how many individual frames 
cover each pixel of the final mosaics. 


2. SAMPLE AND OBSERVATIONS 

To assemble the S^G sample we made use of the 
HyperLEDA database (jPatiirel et al.l 120031 1. We se¬ 
lected all galaxies with radio-derived radial velocity 
V < 3000 km/s, which translates into a distance cut 
of d ;< 40 Mpc for Hq = 71kms“^ Mpc“^. We only 
considered galaxies with total corrected blue magnitude 
"^Scorr < 15.5, blue light isophotal diameter D 25 > I.O', 
and Galactic latitude |&| > 30°. The final sample af¬ 
ter applying all these selection criteria comprises 2352 
galaxies. Using radio-based velocities biases the sample 


4. PIPELINE 2: OBJECT MASKS 

For each galaxy in the sample, P2 creates masks of con¬ 
taminant sources like background galaxies, foreground 
stars and artifacts. Such masks are necessary for sub¬ 
sequent analysis such as background measurements and 
surface photometry (P3), 2D image decompositions (P4) 
and mass map generation (P5). 

The mask creation procedure consists of two steps: 
(i) automatic generation of initial masks, and (ii) visual 
check and editing to create final clean masks. A first 
set of masks is automatically generated by running SEx- 
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tractor (|Bertin fc Arnoutill996D on the 3.6 and 4.5/xm 
images; the resulting segmentation maps constitute our 
initial raw masks. Masking sources is relatively easy in 
those areas of an image with little or no emission from the 
main target galaxy. However, it becomes more compli¬ 
cated on top of the galaxy, where we need to avoid mask¬ 
ing regions belonging to the galaxy itself. Therefore, for 
each galaxy and band we create three automatic masks 
with high, medium, and low detection thresholds in SEx- 
tractor that control how aggressively different sources are 
masked on the main body of the galaxy. 

During the visual quality check step, for each galaxy 
we first choose the best mask among the three ones with 
different thresholds. This best mask is then visually in¬ 
spected and edited by hand, masking additional sources 
missed by SExtractor and unmasking any regions of the 
galaxy that may have been masked by mistake. In partic¬ 
ular, we often have to manually add to the masks the ex¬ 
tended haloes and diffraction spikes of bright stars. Con¬ 
versely, certain sources on the target galaxies like bright 
clumps or the ansae of bars are sometimes picked up by 
SExtractor, and need to be excluded from the masks by 
hand. This editin g process is perfo rmed with a custom 
code described in ISalo et al.l (1201511 . This IDL routine 
displays side by side the original and masked images, 
and provides several geometric shapes that the user can 
employ to interactively mask or unmask regions. The 
editing is first manually done on the 3.6 /rm images, and 
then the editing is automatically transferred to the 4.5/im 
ones, checking that artifacts remain properly masked. 
Thus, we end up with two final, edited masks for each 
galaxy, one for each band. A sample mask is shown in 
Fig. [11 

The masking process is to some extent subjective, es¬ 
pecially for the faintest sources, due to the lack of color 
information. Theses masks are mainly intended to re¬ 
move extraneous sources that would significantly con¬ 
taminate and distort our radial profiles and integrated 
magnitudes. Therefore, users interested in specific struc¬ 
tures such as very faint HII regions or globular clusters 
around our galaxies should not blindly use our masks 
without double-checking the nature of these sources with 
ancillary multi-wavelength data. 

5. PIPELINE 3: SURFACE PHOTOMETRY 

In this section we describe in detail the inner work¬ 
ings of Pipeline 3. Briefly speaking, P3 first determines 
the background level of the images. It then locates the 
centroid of the galaxy and performs different sets of el¬ 
lipse fits. Finally, it measures several properties such as 
asymptotic magnitudes, isophotal sizes and concentra¬ 
tion indices. 

5.1. Sky measurement 

A careful determination of the background level and 
noise is essential to perform reliable surface photome¬ 
try. This is often done by placing several “sky boxes” 
around each galaxy and measuring the corresponding 
background statistics in them. With large samples such 
as S'^G it is desirable to implement this process in a way 
that is automatic yet flexible. 

Pipeline 3 automatically places several sky boxes 
around each galaxy (see Fig.|T]). This is done by defining 



Fig. 1. — Sample object mask (red) and sky boxes (black) overlaid 
on the 3.6 pm image of NGC0936. Note that the sky boxes are 
automatically grown in the radial direction to ensure that they all 
contain the same amount of unmasked pixels. 


two concentric and adjacent elliptical annuli that sur¬ 
round the entire galaxy. Each ring is then azimuthally 
subdivided into 45 sectors or boxes. Since these boxes 
will contain in general a certain amount of masked pix¬ 
els, each box is grown radially outwards until each one 
contains 1000 unmasked pixels. We then measure the 
median sky level and local rms inside each box, as well 
as the large-scale rms between the sky levels of all boxes. 
While the number of 1000 pixels per box is to some extent 
arbitrary, it is chosen to provide a reliable measurement 
of the local rms, while at the same time yielding boxes 
that are small enough to be easily accommodated within 
our images. Note also that before using the P2 masks, 
we first grow the masked areas by 2 pixels to make sure 
that the faint wings of the PSF do not contaminate our 
sky measurements. 

The sky boxes are initially placed by default at 2 x R 25 
from the galaxy center, but this value is modified as 
needed for each galaxy in order to ensure a proper back¬ 
ground subtraction. To do this, we compare the values 
between the inner and outer rings to make sure that there 
are no significant differences, which could be a telltale 
sign of contamination from the galaxy itself. We double 
check this by verifying that the growth curve is flat (see 
Section 1531 ). We also check that the boxes are not too 
close to the frame edges (which are noisier as a result 
of the dithering pattern), or that they do not fall in the 
adjacent frames (which may have a somewhat different 
background value). In cases with a complicated back¬ 
ground structure, the pipeline allows one to manually 
distribute the sky boxes as deemed appropriate. 

The distribution of measured background levels in the 
whole S'*G sample is shown in Fig. |3 The histograms 
peak at /r ~ 24 and 23.5 AB mag arcsec“^ for the 3.6 and 
4.5 /rm bands, respectively. 
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p (AB mag/arcsec^) 

28 27 26 25 24 23 22 



Fig. 2. — Distribution of the sky level within the full S^G sample. 
Red and blue histograms show galaxies from the warm and archival 
(i.e. cryo) groups, respectively. Solid lines correspond to 3.6 pm, 
and dashed ones to 4.5 pm. 

It is illustrative to verify whether our measured val¬ 
ues of the background level agree with theoretical expec¬ 
tations. The background in the S^G images is almost 
entirely dominated by zodiacal light (or “zodi”) coming 
from dust grains in the ecliptic plane. Both the ther¬ 
mal emission of the grains and the scattered sunlight 
contribute to the zodi. At 3.6 and 4.5 pm, both contri¬ 
butions are roughly equal. Thermal emission from in¬ 
terstellar cirrus in the Milky Way amounts to merely 
l%-3% of the background in our images, as it peaks at 
much longer wavelengths (and the S'^G sample specifi¬ 
cally avoids the Galactic plane anyway). Contamination 
by the unresolved Cosmic Infrared Backgr ound is negli¬ 
gible at these bands compared to the zodi (|Hauser et al.l 
[l99p). 

Figure |3] shows how the observed sky level varies with 
the ecliptic latitude of each galaxy. As expected, since 
the primary source of background emission at these wave¬ 
lengths is zodi, the distribution clearly peaks in the eclip¬ 
tic plane. For each galaxy we retrieved from the FITS 
header the background level predicted for its Galactic 
coordinates and epoch of observation0E3 The green 
curves in Fig. [3] are the upper and lower envelopes of 
these predicted values. We can see that both the depen¬ 
dency with ecliptic latitude and the range in background 
level at fixed latitude nicely agree with our measured 
values. 

Interestingly, about 30 WARM galaxies exhibit back- 

For details on the Spitzer background estimator, see 
http://ssc.spitzer.caltech.edu/warmmission/propkit/som/bg/ 
The predicted zodi at the date and coordinates of observa¬ 
tion is stored in the FITS header keyword ZDDY_EST. Part of this 
zodi is removed when the Spitzer automatic pipeline subtracts a 
skydark frame; the estimated amount of removed zodi is given by 
SKYDRKZB. The final level of zodi that remains in our images is 
therefore ZQDY_EST - SKYDRKZB. 



-80-60-40 -20 0 20 40 60 80 -80-60 -40-20 0 20 40 60 80 


bed (beg) 


bed (beg) 


Fig. 3.— Background level measured in our fully-reduced images 
as a function of the ecliptic latitude of each galaxy. Red and blue 
dots correspond to warm and archival galaxies. The green curves 
delimit the distribution of the zodi brightness predicted by the 
Spitzer background estimator. 



Fig. 4. — Two sample images exhibiting a diffuse background 
artifact: NGC5597 (left) and UGC08507 (right). The smooth pat¬ 
tern is seen both in the main frame (top) and in the flanking one 
(bottom). A histogram equalization has been applied to exagger¬ 
ate the brightness of this artifact. This structure only appears in 
~ 30 of our galaxies. 


ground levels at 3.6 /tm well above the predicted values. 
A visual inspection of these images revealed a recurrent 
diffuse artifact, both in the main frame where the galaxy 
is and in the flanking one (Fig. 2]). These few galaxies 
were observed in a narrow time window during August- 
September 2009, very shortly after the Spitzer’s warm 
phase began. This smooth background artifact may be 
therefore due to the detector temperature and bias not 
being settled yet at that time. 

Besides the background level itself, it is also impor¬ 
tant to characterize the background noise at different 
spatial scales, as this determines our ability to detect and 
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Limit n (AB mag/arcsec^) 

27 26 25 24 23 icr 

26 25 24 23 22 5o 



Limit ji (AB mag/arcsec^) 

29 28 27 26 25 lo 

28 27 26 25 24 5a 



log (large scale rms) (MJy/sr) 


Fig. 5.— Distribution of local, pixel-to-pixel noise (left) and large 
scale background noise (right) in the S^G images. WARM and 
ARCHIVAL galaxies are shown in red and blue, respectively. Solid 
and dashed lines correspond to 3.6 and 4.5 fim. The corresponding 
surface brightness limit at 5 (j and Icr is shown at the top. 

measure faint structures at the outskirts of our galaxies. 
Figure [5] shows the distribution of the local, pixel-to-pixel 
noise in our images. The histogram for the warm galaxies 
is much narrower than that for the archival ones, which is 
expected given that the warm galaxies were imaged with 
the same observing strategy, whereas the archival come 
from a variety of different programs. In general we reach 
a local surface brightness sensitivity per pixel of ~ 24 
and ^ 26 AB magarcsec”^ at 5cr and Icr, respectively. 

Note, however, that these are values for individual pix¬ 
els. When measuring the flux of a given extended source, 
this noise component scales down with the square root 
of the number of pixels in the region where the mea¬ 
surement is being performed. In particular, when mea¬ 
suring radial profiles at the outskirts of galaxies, this 
local noise component becomes negligible compared to 
the large scale background noise, which is given by the 
rms between the median sky values measured in the 
different boxes. As shown in Fig. [Sj this large scale 
component peaks at 26-26.5 AB mag arcsec“^ at a 5cr 
level, and ^ 28 ABmagarcsec“^ at Icr. As a reference, 
/is.6 = 27 AB magarcsec”^ corresponds to a stellar mass 
surface density of I Mr?, adop ting = 0.53 

as me asured bv lEskew et all (1201211 : see also lMeidt et al.l 

(l2?fn . 


5.2. Radial profiles 

Once the background level and noise have been mea¬ 
sured, our pipeline proceeds to perform surface photom¬ 
etry on the images. We first use the IRAlIU task IM- 
CENTROID to find accurate coordinates for the center of 
the galaxy. We the n run the task ellipse (j.Iedrzeiewskil 
Il987t iBnskol 1199611 to obtain radial profiles of surface 
brightness {p), ellipticity (e) and position angle (PA). 
We perform three separate runs of ellipse with differ¬ 
ent settings: 

1. Fixed center; free e and PA; radial resolution 

Ar = 6". We use these fits with a coarse radial in¬ 
crement to derive e and PA in the outer parts of the 
galaxy, thus defining a global shape and orientation 
for each object. For each band, we provide val¬ 
ues of e and PA at two levels of surface brightness, 

IRAF is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 



Fig. 6 .— Sample ellipse fit results for NGC0936 with free e and 
PA, and a radial increment Ar = 2”. The images to the left 
correspond to the 3.6 pm band. In the bottom one we have over¬ 
laid the ellipses fitted by our pipeline. The red solid ellipse cor¬ 
responds to p 3 ,g = 25.5 AB mag arcsec“^, and the dashed one to 

26.5 AB mag arcsec“^. The blue ellipse is the optical D 25 one from 
HyperLeda. In the radial profiles, offsets have been applied to the 

4.5 pm data as indicated in the legend, to avoid overlapping. The 
right vertical scale in the top panel marks the stellar mass surface 
density for the 3.6 pm profile. 


pL = 25.5 and 26.5 AB mag arcsec“^. After testing 
how sensitive these values are to variations in the 
sky subtraction, input fitting parameters, amount 
of masked objects, etc., we recommend using the 
25.5 values as they are more stable. 

2. Fixed center; e and PA fixed to the values at 
25.5 AB magarcsec”^; Ar = 2". These fits have 
a finer radial resolution that matches the IRAC 
PSF at these wavelengths. By keeping e and PA 
fixed and equal to the global outer values, these 
fixed-fits are ideal to measure disk scale-lengths, 
disk break radii, to perform ID bulge-disk decom¬ 
positions, etc. We also employ these profiles to 
measure the integrated magnitude of each galaxy 
from the growth curve (Sect. 15.41) . 

3. Fixed center; free e and PA; Ar = 2". These 
free-fits with a fine resolution are well suited to 
study in detail the structural properties of features 
such as bars, which leave very characteristic signa¬ 
tures in the radial profiles of e and PA. An example 
of these fits is shown in Fig. [51 

For simplicity and to improve the robustness of the fits, 
once we measure the central coordinates of galaxy, these 
are kept fixed during the ellipse fitting. However, the 
pipeline allows one to leave the center as a free param¬ 
eter too. We have in fact performed such fits to study 
individual galaxies with offset bars (e.g. NGC3906, de 
Swardt et al. 2014, in prep.). 

We normally begin the ellipse fitting at an intermediate 
radius, typically between 0.5 — 1 x i? 25 , using the optical e 
and PA from HyperLEDA as input guesses for the fitting 
routine. Starting at that initial radius, the ellipse fitting 
proceeds first outwards and then inwards. 

It should be noted that the ellipse algorithm was 
specifically designed for galaxies with a smooth radial 
brightness distribution. It therefore works correctly for 
elliptical galaxies, but it does not cope well with fea¬ 
tures such as spiral arms or rings that can significantly 
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modify the local luminosity gradient. As a result, el¬ 
lipse can sometimes stop fitting at a given radius in 
spiral galaxies. Tuning the input parameters can some¬ 
times fix this problem and provide a good fit, but in 
general this is a blind trial-and-error solution with an 
unpredictable outcome. Some authors have in fact cir¬ 
cumvented this issue by automatically running ellipse 
tens or even hundreds of times on each image, varying 
the initial fi tting parameters each time until a good fit 
is obtained (IJogee et al.ll2004ll . While this approach can 
be practical in small images of distant galaxies, it is pro¬ 
hibitively time-consuming for the large images in S^G, 
where a single fit can take up to several minutes (and 
this is further complicated by the large number of galax¬ 
ies in our sample). We thus decided to implement an 
optimized procedure that is better tailored to our needs. 
After each run of ellipse, the pipeline looks for any ra¬ 
dial interval where the fit did not converge, and it then 
refits that particular radial interval alone. The process 
is then iterated several times to ensure that no radial 
gap is left unfit. In general most of our galaxies can be 
properly fit with a single ellipse run, but this automatic 
iterative process was of great help to handle troublesome 
cases that would have otherwise required considerable 
manual work. 

In any case, for each galaxy we always visually inspect 
the output of the pipeline: we overlay the fitted ellipses 
on the galaxy image, and we plot the radial profiles of /r, 
e and PA, making sure that the results accurately probe 
the different structures within the galaxy. 

Dwarf galaxies can be particularly challenging for our 
pipeline. They are often partially resolved into stars in 
our images, and their clumpy and patchy appearance can 
sometimes fool the ellipse fitting algorithm. For these 
galaxies we recommend not to overinterpret any signa¬ 
ture in the e and PA profiles without inspecting the im¬ 
ages themselves. Since these galaxies normally lack well- 
defined large-scale structures such as bars, the profiles 
with fixed e and PA are more suitable in these cases. 

We have released all these ellipse fits in the form of 
ASCII tables containing the full output of the ellipse 
task. These tables not only include radial profiles of fi, e 
and PA, but also other quantities such as the harmonic 
deviations from a perfect ellipse, which are commonly 
used to quantify the boxiness/diskiness of the isophotes 
at diff erent radii (e.g. ICarteil[l978t iKormendv fc Bended 
[ 19 ^ . 

Apart from the usual output from ellipse, we also in¬ 
clude in our ASCII profiles additional columns that are 
specific to the IRAC data used here. In particular, the 
IRAC photometry needs to be corrected for the extended 
wings of the PSF and the diffuse light that is scattered 
throughout the detector. Here we rely on the extended 
source aperture correction provided in the IRAC Instru¬ 
ment HandboolF^. Given an elliptical aperture with ma¬ 
jor and minor radii a and 5, if Fhbs is the total observed 
flux inside such an aperture, the corrected flux is given 
by: 

7^corr(req) = Aobs(Peq) X -b C) (1) 

In this expression req = is the equivalent radius of 
the elliptical aperture, in arcseconds. The coefficients A, 


B and C are equal to 0.82, 0.370 and 0.910, respectively, 
at 3.6/rm, and 1.16, 0.433 and 0.94 at 4.5/xm. 

Similarly, if Jobs is the surface brightness along an 
isophote (rather than the total flux inside that radius 
as before), the aperture-corrected surface brightness can 
be obtained by performing a series expansion on the pre¬ 
vious equation: 

7corr(l"eq) —7obs(l"eq) X (Ac C) 

- ABr^“^e“’'^^'>Fobs(i’eq)/(27r) (2) 

Our ASCII tables include radial profiles with and with¬ 
out these aperture corrections. These corrections are es¬ 
timated to be uncertain at a 5%-10% level. We have 
always applied the extended source aperture correction 
at all radii in our profiles. However, the reader should 
keep in mind that the point source corrections might be 
more suitable at very small radii (r < 8-9"), especially if 
the nucleus is bright and compact. 

Finally, all the released profiles have been corrected 
for foreground ex t inctio n, using the color excess map 
of iSchleeel et al.l (I1998D a nd th e Milky Way extinc- 
tion curve of ILi fc Draini (1200111 . Nevertheless, at 3.6 
and 4.5/im this extinction correction is typically around 
0.005 mag. 

5.3. Error analysis 

The pipeline also estimates the uncertainty in fj,, e and 
PA at different radii. For e and PA we rely on the er¬ 
rors determined by the ellipse task, which result from 
the internal errors in the harmonic fit (lBuskolll996ll . As 
for /i, the error provided by ellipse is derived from the 
rms of the pixel values along each isophote, so it mostly 
reflects azimuthal variations of the stellar emission at a 
given radius rather than the true uncertainty of the mean 
surface brightness at that radius. We therefore opted to 
measure the error in /x ourselves. 

If / is the incident pixel intensity, prior to sky subtrac¬ 
tion, and /sky is the sky level, then the surface brightness 
is given by: 

/x = -2.51og(/-/sky)+/f (3) 

where K includes the magnitude zero-point plus any 
other global multiplying factor, such as the aper¬ 
ture corrections mention ed before. Followi n g the 
met hodology detailed in Gil de Paz fc Madord (|2005[1 
and IMunoz-Mateos et al.1 ( 2009all . the uncertainty in fj, 
can be estimated as: 

Am= ^(A/F)2+ (A/^ + A/i^) (4) 

The error in K is dominated by the uncertainty in 
the aperture corrections, and is of the order of 5-10% as 
explained in the previous section. Since any change in 
K will merely translate into a global offset of the radial 
profiles, the values of A/x in our ASCII tables and plots 
do not explicitly include the contribution of AK. 

We computed A/, which is the error in the incident 
pixel intensity, by assuming poissonian statistics: 


' http://irsa.ipac.calt6ch.edu/data/SPITZER/docs/irac/iracinstruinenthaiidbook/30/li/ = 


Heff A^isophote 


( 5 ) 
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where ^eff is the effective gain and -/Visophote is the number 
of pixels within each isophote. The effective gain, in 
turn, was derived from the nominal detector gain and 
the exposure time of each pixel according to the weight 
maps. 

The term A/gky represents the uncertainty in the back¬ 
ground level, and was obtained as: 


A/sky — 


^ local 


iVi 


isophote 


+ max 1 - 


local 


iVbo 


( 6 ) 


where uiocai and criarge are the local and large scale back¬ 
ground rms, respectively, as described in Section 15.11 
Here Aisophote is, again, the number of pixels along a 
given isophote, and Abox is the number of pixels inside 
the boxes used to measure the background level. 

The first term in Eq. [S] reflects the contribution of 
the local, pixel-to-pixel noise to the final error in the 
surface brightness. This contribution scales down as 
1/y^Aisophote, so it becomes negligible at large radii, 
where the flux is averaged over a large number of pix¬ 
els. It is the second term, the large-scale noise, the one 
that dominates the error budget at large radii. Note 
that the measured rms between the sky values in the dif¬ 
ferent boxes is partly contaminated by the local noise, 
in the sense that even if there were no true large-scale 
fluctuat ions, the rms between the sky boxes would be 
ciocai/VNbox, on average. The correction term in Eq. |6] 
accounts for this. 


5.4. Asymptotic magnitudes 

From the surface photometry we measure the asymp¬ 
totic magnitude of each galaxy by extrapolating the 
growth curve to infinity. To do this we use the profiles 
with fixed e and PA and 2" resolution. As an example, 
the bottom panel of Figure [3 shows the 3.6/rm growth 
curve of NGC0936, which rises fast at the center and 
then slowly approaches a flat regime. To determine this 
asymptotic value, we first compute the gradient of the 
growth curve at all radii. In the outer parts of the galaxy, 
the gradient at a given radius and the magnitude inside 
that radius follow a linear trend (top panel). We then 
apply a linear fit whose y-intercept (the magnitude for a 
flat gradient) is by definition the asymptotic magnitude. 

Note that growth curves are often used in shallow im¬ 
ages to recover the flux buried beneath the noise in the 
outskirts of galaxies. This is not the case in the S'*G im¬ 
ages, which are deep enough to always comfortably reach 
the flat regime of the growth curve. Therefore, in prac¬ 
tice there is no extrapolation involved when fitting the 
growth curves of our galaxies. In this regard, the error 
bugdet in the surface photometry is dominated by the 
large-scale rms, and not by any possible non-zero slope 
in the growth curve. 

Once we have determined the asymptotic apparent 
magnitudes, we derive the absolute ones using redshift- 
independent distances from NED whenever available 
(79% of the sample), and redshift-dependent ones other¬ 
wise, assuming Hq = 71kms“^ Mpc“^. For those galax¬ 
ies with more than one redshift-independent distance 
measurement (61% of the whole sample), the typical 
rms between the different measured distances is ~ 15%, 
which translates into an error of ^ 0.3 mag in the ab¬ 
solute magnitude. While this is not strictly speaking a 



Fig. 7. — Illustration of how the asymptotic magnitudes are 
computed. The main panel shows the 3.6 /rm growth curve of 
NGC0936, that is, the integrated magnitude as a function of ra¬ 
dius. The small inset shows the magnitude as a function of the 
magnitude gradient for the data-points in the flat portion of the 
growth curve. The red line represents a linear fit to both quantities. 
The asymptotic magnitude is the y-intercept of this fit, namely, the 
magnitude for a flat gradient (blue circle). 


true formal uncertainty, for such nearby galaxies it pro¬ 
vides a more realistic estimate of the actual distance error 
than the uncertainty in the radial velocity, which is just 
^ 0.3% for our galaxies. 

In Fig.|8]we plot the 3.6^m absolute magnitude for all 
our galaxies as a function of their optical morphological 
type as given by HyperLEDA (the 4.5/rm values follow 
a very similar trend). The conversion to stellar mass 
in the ri ghtmost vert i cal ax is was done via the M*/L 
value of feskew et al.l (I2012D . At the late-type end of 
the sequence, the stellar mass rises monotonically from 
^ 5 X 10® Mq for irregular galaxies to 1 — 2 x 10^° Mq 
for Sc ones. The trend then remains remarkably flat out 
to SO galaxies, beyond which the stellar mass increases 
again for elliptical galaxies. The rms in stellar mass for 
each bin of Hubble type is typically between 0.5-0.6dex, 
or a factor of 3-4. 


5.5. Concentration indices 

Gentral light concentration is one of the parameters 
that varies most prominently along the Hubble sequence, 
and was in fact used early to modi fy Hubble’s original 
classification scheme (jMorganI 119581) . Light concentra¬ 
tion and other non-parametric estimators such as the 
asymmetry, dumpiness, second-order moment or Gini 
coefficient, haye been extensiyely used oyer the years to 
qua ntify the morphology of nearby and distant galax- 

ies (iBershady et al.ll200(]l:lConselice^200^lConselice et ahl 


20031: lAbraham et al. Il996bllal. ^03t Lotz et al. 200 


TayionAfaggret^_ 2()07t IMufioz- Mateos et al.1 2009 
Holwerda et al.1 1201 lib In particular. iHolwerda et al l 
(j2013l ) presented a detailed study of the quantitatiye 
morphology of the S'^G galaxies through a suite of com¬ 
monly used non-parametric estimators. 

Here we present the concentration indices deriyed from 
the growth curves measured by our pipeline. Among the 
different definitions existing i n the liter ature, we se ttled 
on Gsi (jde Vaucouleursl[l977l) and Cs 2 (|Kentlll9^ . de- 
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Fig. 8. — Stellar mass and absolute magnitude at 3.6/i.m of all 
galaxies as a function of their morphological type, together 
with the corresponding histograms of both quantities. The solid 
and dashed lines show the median trend and rms in bins of AT = 1. 

fined as follows: 

c,. = ni (T) 

r25 

C'82=51og (8) 

\r 20 J 

where is the semi-major axis of the ellipse enclosing 
x% of the total luminosity of the galaxy. There is some 
discrepancy in the literature regarding how the total lu¬ 
minosity is measured when computing concentration in¬ 
dices; here we simply use the asymptotic magnitudes de¬ 
scribed before. Also, unlike concentration indices derived 
from Sersic fitting, our values are truly non-parametric, 
in the sense that we do not assume any functional form 
for the radial profiles or the growth curves. 

For each galaxy and band we computed both concen¬ 
tration indices. The results are presented in Sec. 16.11 
where we investigate why galaxies with the same stel¬ 
lar mass can exhibit vastly different morphologies and 
concentration indices. 

5.6. Galaxy sizes and shapes 

From the growth curves we have also obtained effec¬ 
tive radii (res) containing half of the total asymptotic 
luminosity, independently at both IRAC bands. Ra¬ 
dial gradients in age, metallicity and internal extinction 
have a negligible impact at these wavelengths; therefore, 
our half-light radii can be safely interpreted as half-mass 
radii. Given that we measured reS along the major axis 
of elliptical isophotes, our values do not need to be cor¬ 
rected for inclination, as is the case in studies employing 
circular apertures. Also, since we do not perform any 
sort of Sersic or similar fitting, the values that we pro¬ 
vide constitute the true effective radius of each galaxy, 
and not the effective radius of a simple model fitted to 
that galaxy. In Sect. 16.21 we discuss in detail the stellar 
mass-size relation of the S^G galaxies for different mor¬ 
phologies. 


The effective radii depend on the level of central con¬ 
centration, in the sense that highly-concentrated galaxies 
tend to have small values of Teff compared to the overall 
extent of those galaxies. Therefore, it is also desirable 
to measure the global sizes and shapes of our galaxies, 
including their outermost parts. Here we follow an RC3- 
like approach and measure the radius, ellipticity and po¬ 
sition angle at 25.5 and 26.5 AB mag arcsec“^, both at 3.6 
and 4.5/rm (that is, four sets of measurements for each 
galaxy). We recommend using the 25.5 values at 3.6^m 
to characterize the overall extent of the S'^G galaxies. 
Even though there is always plenty of detectable emis¬ 
sion beyond that surface brightness level, ellipse fitting 
becomes less reliable and more sensitive to convergence 
issues, object masking and background subtraction. The 
distribution of isophotal sizes as a function of stellar mass 
and morphology is presented in Sect. 16.21 

Regarding the outer isophotal ellipticities, they can be 
used as a photometric proxy for the disk inclination in the 
case of spiral galaxies. We have in fact used the P3 ellip¬ 
ticities to deproject our galaxies and rec over the intrin¬ 
sic shapes of rings and bars, for i nstance (|Comer6n et al.l 
120131 : iMuhoz-Mateos et al.ll2013ll . 

Note, however, that this is complicated by the fact that 
(a) disks may not be perfectly circular when viewed face- 
on, and (b) the intrinsic vertical thickness of disks and 
bulges/haloes will dist ort the outermost is ophotes at high 
inclinations (see, e.g.. lRvd^l2004 120061 and references 
therein for a discussion on both issues). The influence 
of the vertical thickness is particularly obvious in our 
deep IRAC images, given that thick disks and spheroids 
are primarily composed of old stellar populations, which 
show up prominently at these wavelengths. Users should 
therefore keep this in mi nd when studying high ly-inclined 
galaxies in our sample (iGomeron et al.lT2011ll . We will 
present inclinations derived from multicomponent 2D fit¬ 
ting of the galaxy images in a forthcoming paper (Salo 
et al. 2014, in prep.). 

6 . RESULTS AND DISCUSSION 

6.1. Variations in concentration and moirphology at 
fixed stellar mass 

Figure [5] shows the distribution of the Cs 2 index at 
3.6/j,m as a function of the stellar mass for all S'^G galax¬ 
ies. FigureHnidisplays the same information in histogram 
form. None of these plots change noticeably when using 
C 31 and/or the 4.5/im data instead. 

The concentration index distribution exhibits a pro¬ 
nounced narrow peak centered at Cs 2 2.8-3, and a 
much broader one peaking at Cs 2 ^ 4.5-5. The first 
group corresponds to disk-dominated galaxies; indeed, 
the theoretical concentration index for a perfectly expo¬ 
nential profile is 2.8. The second group comprises bulge- 
dominated galaxies and ellipticals. The “L”-shaped dis¬ 
tribution of concentration vs. stellar mass shown in 
Fig. [9] nicely describes how the spatial distribution of old 
stars within galaxies varies across the Hubble sequence. 
For galaxies later than Sc, the concentration index re¬ 
mains roughly constant at the value expected for disk- 
dominated galaxies, with some scatter but no obvious 
dependence on the stellar mass. But for galaxies earlier 
than Sc the situation changes, and now the concentration 
index ranges anywhere from a disk-dominated profile to 
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a highly concentrated on e . A similar distribution was 
found bv iScodeggio et ahl (I2002D in their analysis of the 
iJ-band photometry of galaxies in nearby clusters. They 
found that galaxies with Lh < exhibit a small 

and constant concentration index, as well as blue optical- 
IR colors. At Lh > galaxies were found to span 

a large range of concentration indices, with redder colors. 

To better quantify the structural variety of the S^G 
galaxies, in Fig. [TT] we show the median profiles after 
grouping the galaxies in bins of stellar mass and concen¬ 
tration. Based on the fact that Cs 2 = 2.8 for a pure 
exponential profile and C ^2 = 5.3 for a de Vaucouleurs 
one, we use these values to define four bins of mass con¬ 
centration, as shown in Fig. 1111 blue and green curves are 
the median profiles of disk-dominated galaxies, with the 
blue ones being less concentrated than the green ones. 
On the other hand, orange and red curves correspond to 
bulge-dominated galaxies, with the latter being the most 
concentrated ones. 

In galaxies more massive than 10^^ Mq, the observed 
differences in concentration at a fixed mass are driven by 
variations in the global radial stellar structure over scales 
of many kpc, all the way from the center of the galaxies 
to their outer parts. However, for galaxies less massive 
than 10^° Mq the situation is different: differences in 
concentration between the blue and green curves (both 
disk-dominated) are mostly due to the presence or ab¬ 
sence of central bright features smaller than 1 kpc; the 
disks at r > 1 kpc are otherwise very similar in terms of 
surface density and slope. 

In the following two subsections we investigate in more 
detail why galaxies with the same stellar mass and mor¬ 
phological type present such varied radial structure and 
concentration indices. 

6.1.1. Spread in concentration above lO*^® M© 

According to Fig. [9l the structural transition between 
low- and high-concentration galaxies occurs at stellar 
masses between 10^° and 10^^ Mq. Previous studies 
based on optical imaging and spectroscopy have found 
abrupt transitions in the stellar ages and star forma¬ 
tion h i stories in the same mass in terval (|Kauffmann et ^ 
120031: iBrinchmann et id] l2004f l. with stellar popula¬ 
tions suddenly becoming older when going from low- 
mass, disk-dominated galaxies to more massive, highly- 
concentrated ones. 

To quantify changes in the stellar populations of the 
S^G galaxies we made use of GALEX measurements 
in the far and near ultraviolet. Bouquin et al. (2015, 
submitted) measured FUV and NUV surface photom¬ 
etry and asymptotic magnitudes for the S^G galaxies 
following the same methodology described here for the 
Spitzer images. Here we use their measurements to de¬ 
rive extinction-corrected star formation rates (SFR). 

For each galaxy we first estimated the internal dust 
extinction via the integrated FUV —NUV color, follow¬ 
ing the prescription calibrated by iMuhoz-Mateos et al.l 
(|2009bH on nearby galaxies. While this recipe is a good 
proxy for internal extinction in normal star-forming spi¬ 
rals, it overestimates the true extinction in elliptical 
galaxies, where a significant fraction of the observed UV 
reddening is due to the intrinsicially old and red stel¬ 
lar populations. It also overestimates the extinction in 
Sdm-Irr galaxies, which tend to have somewhat redder 


UV colors than other low-mass disks with similar dust 
content, possibly due to differences in the extinction 
law, dust geometry, or a bursty star formation activ¬ 
ity. Therefore, we only relied on the FUV—NUV color to 
determine the internal extinction i n the disk-dominated 
galaxie s of our sample. Following iMnhoz-Mateos et al.l 
(l2009bli , for elliptical galaxies we adopted a constant ex¬ 
tinction in the FUV of 2 mags, and 0.5 mags in the case 
of Sdm-Irr galaxies. The extinction-corrected FUV lumi- 
n osities were then c onverted into SFR via the calibration 
of iKennicnttI (jl998l l. Note that in elliptical galaxies most 
of the UV light is not associated to recent star formation, 
so t heir true SFR is lower than the value resulting from 
the iKennicutti (I1998D prescription. 

The data-points in Fig. are colored according to their 
specific SFR (sSFR, the SFR per unit of stellar mass). 
The transition from low- to high-concentration galaxies 
in the mass regime between 10^° and 10^^ Mq is ac¬ 
companied by a decrease of one order of magnitude in 
the sSFR, which drops from ^ 10“^°'^yr“^ in massive, 
low-concentration galaxies to ^ 10“^^'^yr“^ in high- 
concentration galaxies with the same stellar mass. The 
same decrease in sSFR is found in spectroscopic studies, 
where the sSFR is estimated from spectral measurements 
such as the break at 4000A or the Ha equivalen t width 
(|Kauffmann et al]l2003t IBrinchmann et al1l2004|) . 

This drop in sSFR with increasing concentration does 
not occur all at once for all morphological types, though. 
Figure[H]demonstrates that it is only in SO/a galaxies and 
earlier types that the sSFR drops below ^ 10“^^yr“^ in 
highly-concentrated objects. In these galaxies the high- 
concentration and red colors result from the dominant 
contribution of large bulges and stellar haloes to the total 
light. But in the morphological range between Sa and 
Sbc, highly-concentrated galaxies still exhibit a relatively 
high sSFR of ~ 10“^°yr“^, similar to the sSFR of less 
concentrated galaxies in the same morphological regime. 

What is the physical nature of these galaxies with a 
concentrated stellar mass distribution but high sSFR? A 
visual inspection of these objects reveals that they are a 
mixed bag of galaxies with very different physical prop¬ 
erties, which can be nevertheless sorted out into three 
broad categories fFig. [T^: 

1. Barred galaxies with prominent nuclear rings, such 
as NGC 7552, NGG 4593 or NGG 1365. The 
strong bars in these galaxies drive gas inwards, but 
this gas eventually stalls in a nuclear ring, most 
likely associated with the Inner Lindblad Reso - 
nance (IRegan et al.lflMl IShddi et al.ll2nnnl [209^ . 
This shrinks r 2 o, the radius of the isophote contain¬ 
ing 20% of the total luminosity of the galaxy, thus 
increasing the concentration index Cs 2 in these 
galaxies. Note that given the intense star forma¬ 
tion activity in nuclear rings and their elevated 
gas and dust content, non-stellar sources can lo¬ 
cally dominate the 3.6 and 4.5 /im emission in these 
rings. Nevertheless, the clean stellar mass maps 
resulting from our P ipeline 5 (jMeidt et al.l 120121 : 
iQnereieta et al.l[2014f ) still reveal a prominent cen¬ 
tral concentration of stellar mass in these galax¬ 
ies, indicating that bar-driven gas inflows play a 
key role in the secular assembly of (pseudo- ) bulg es 
(jkormendv fc Kennicutdl2004 iSheth et al.ll2005r) . 
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Fig. 9. — Concentration index Cs 2 at 3.6/rm as a function of total stellar mass. The gray datapoints show the full S'^G sample, whereas 
the colored ones correspond to galaxies with different Hubble types, as indicated in each panel. The color indicates the extinction-corrected 
specific SFR, computed from GALEX and Spitzer data as explained in the text. As a reference, in the top-left panel we plot with dashed 
horizontal lines the theoretical concentration indices for Sersic profiles with different indices n. 


2. Interacting systems such as NGC 2782, NGC 7714 
or NGC 5534. These galaxies often exhibit up- 
bending radial profiles, with a steep inner disk fol¬ 
lowed by a flatter outer one, normally accompanied 
by tidal features in the outermost regions. This 
decreases r 2 o and increases rgo, yielding high con¬ 
centration indices. Numeric al simulations of mi - 
nor mergers such as those bv lYounger et all (120071 1 
can reproduce this radial structure: the interac¬ 
tion leads to gas inflows that contract the inner 
radial profile, while at the same time the outwards 
transfer of angular momentum expands the outer 
disk. Also, the sSFR is temporarily enhanced in 
these galaxies during the interaction, with promi¬ 
nent star forming knots that are clearly visible in 
the GALEX images. 

3. Galaxies with compact bulges and smooth ex¬ 
tended disks, often undisturbed as in NGC 3642. 
This again leads to a small r 2 o and large rso- 
While these extended outer disks often have low 
surface brightness in the infrared, they can be 


much brighter at UV wavelengths, as happens with 
NGC 3642. These so-called extended UV disks are 
typically found in ^ 10% - 20% of neyby galaxie s 
(IGil de 'Paz et al.l[20M iThilker [2001 . 

but the origin of these galaxy-wide star formation 
eyents is still under debate. 

In brief, for stellar masses aboye 10 ^'’Mq, galaxies 
with a centrally concentrated radial distribution of old 
stars are either (a): quiescent galaxies with a prominent 
stellar bulge/halo, or (b) star-forming galaxies with cen¬ 
tral mass concentrations due to a yariety of mechanisms 
such as bar- or merger-driyen inflows. 

6.1.2. Spread in concentration below M© 

Figure [9] shows that galaxies with stellar masses below 
10^° Mq still exhibit a considerable scatter in the concen¬ 
tration index at any giyen mass. While all these galax¬ 
ies haye concentration indices broadly consistent with a 
disk-dominated prohle, the actual yalues can range any¬ 
where from Cs 2 ~ 2 to ^ 3.5. Moreoyer, this scatter 
in Cs 2 does not seem to be correlated with changes in 
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Fig. 10. — Histograms of the concentration index Cs 2 at 3.6/i.m 
for the S^G galaxies. The black histogram is the same in all panels, 
and corresponds to the full sample. The red and blue histograms 
show galaxies grouped in bins of morphological type and log M*, 
respectively. 
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Fig. 11. — Median surface brightness profiles of the galaxies 
in bins of log Mq , as indicated by the numbers at the top-right of 
each panel, and concentration index C% 2 ', according to the color 
scheme in the legend. Prior to computing the median, the profiles 
were corrected for inclination and resampled to a commom physical 
scale in kpc. At least five galaxies were required in each bin to 
derive the median. The errorbars in each panel show the ±1 (j 
deviation around the median profile in each case. 

the sSFR, as was the case in the more massive galaxies 
analyzed before. 

To ascertain the origin of these structural variations in 
low mass disks, we visually inspected galaxies with the 
same stellar mass but extreme values of concentration. 
Figure [13] displays several representative examples. 

The top two rows show very low-mass dwarf disks, all 
of them with stellar masses between 10® — 10®'® Mq. The 
objects in the first row are characterized by very low 
concentration indices of Cs 2 ~ 2.0 — 2.2 (smaller than 



Fig. 12.— Three sample S'^G galaxies with high concentration 
indices at 3.6 fim {Cs 2 ~ 5.0 — 5.5) and also high specific SFR (< 
yr”^): NGC 7552, a barred spiral, NGC 2782, an interacting 
disk, and NGG 3642, a galaxy with a compact nucleus/bulge and 
a diffuse extended disk. The top panels show the S'*G images at 
3.6/rm, whereas the b ottom panels display fa lse color FUV-I-NUV 
images from GALEX IIGil de Paz et aT]|2007ll . 



Fig. 13. — First row: galaxies with stellar masses between 
10® — 10®'® Mq and low concentration indices (C 82 ~ 2.0 — 2.2). 
Second row: galaxies with the same stellar mass as those in the 
first row, but higher concentration values (Cg 2 ~ 3.3). Third row: 
more massive disks with masses around 10®'® — 10^® Mq, but low 
concentration. Fourth row: same stellar mass as in the third row, 
but higher concentration. 

for an exponential profile), whereas those in the second 
row all have Cs 2 ^ 3.3 (larger than for an exponential). 
This higher concentration is due to central aggregations 
of stars (usually resolved in our images) that are ab¬ 
sent in the first group of galaxies. GALEX images of 
some of these concentrated dwarfs, like NGC 0059 and 
NGC 3738, reveal bluer UV colors in the central parts 
than in the outskirts, possibly hinting to an outside-in 
formation for these objects. 

The third and fourth rows in Fig. |T3| correspond to 
more massive disks, with stellar masses between 10®'® — 
10^® Mq. Galaxies in the third row present Cs 2 ^ 
2.0 — 2.2, while those in the fourth row have Cs 2 "^3.3. 
In this case the higher concentration results from bright 










































Surface photometry of galaxies 


13 


bars as in NGC 0600 or NGC 6140. Unbarred galaxies 
can be also concentrated if they host bright and com¬ 
pact (pseudo-)bulges, as in NGC 3344. In contrast, low 
concentrated disks in this mass range either have weakly- 
defined bars (NGC 5147) or no bars at all (NGC 5949). 
Moreover, some galaxies such as NGC 4353 have bright 
inner disks that lead to central plateaus in their radial 
profiles, leading to concentration indices lower than those 
expected for an exponential profile. 

We also find a wealth of asymmetric galaxies which 
may have recently undergone gravitational interactions. 
In some of them, like NGC 2793, the stellar mass is so 
dislodged that the radial surface density is almost flat, 
leading to very low concentration indices. On the other 
hand, galaxies like NGC 1637 exhibit a much milder 
asy mmetry, and therefor e haye higher concentration yal- 
ues (|Zaritsky et al.l[20T^ and references therein). 

In summary, eyen though galaxies with stellar masses 
below 10^° Mq tend to haye low concentration, at any 
giyen mass variations in concentration are still present, 
due to the presence of central star clusters, bars, pseudo¬ 
bulges, or lack thereof. 



56789 10 11 56789 10 11 56789 10 11 


log M, (Mg) 


Fig. 14.— Semi-major radius at /rs.e = 25.5 AB mag arcsec“^ as 
a function of stellar mass of all S^G galaxies (gray dots), grouped 
by their morphological type (black dots). 


6.2. The stellar mass-size relation 

The mass-size relation provides very stringent con¬ 
straints on models of galaxy formation and evolution. 
From a cosmological point of view, the trend between in¬ 
creasing galaxy size and mass reflects, to first order, the 
physical connection between the mass of the dark mat¬ 
ter halo and its angular momentum. Indeed, in a ACDM 
cosmology the characteristic size of a disk should scale as 
(iFall fc EfstathioulligSOt iMo et al1ll998D . In prac¬ 
tice, though, the stellar mass-size relation also depends 
on whether the halo gas conserves its angular momen¬ 
tum as it settles onto the disk, how efficiently that gas is 
later converted into stars, the impact of AGN and stellar 
feedback, and whether internal and external processes 
rearrange stars within the galaxy. 

The mass-size (or luminosity-size) relation has 
been extens ively explored both in the local u n iverse 
(see, e.g.. iKauffmann et al.l l2003t iShen et al.l 2003 : 


Courteau^ ajj_ 2007t iFerndndez Lorenzo et al.f 2013: 
Lange et al.ll2014 ) an d at hig h redsh i ft (e.g. Barden et al l 


2005t iTruiillo et al.l 12004 . 120061: iFranx et al.l I200F . 
van der Wei et al MiiT Here we take advantage of the 
little sensitivity of the IRAC bands to dust extinction 
and M/L variations to derive a self-consistent stellar 
mass-size relation in the local universe. 

In Fig. [T3] we plot R 25 . 5 , the semi-major radius at 
/i 3.6 = 25.5 AB mag arcsec“^, as a function of the to¬ 
tal stellar mass for all S'^G galaxies. We find an obvious 
and expected monotonic trend, in the sense that galax¬ 
ies with higher stellar mass are also larger. However, the 
slope of this trend is not constant, as galaxies of different 
Hubble types populate different areas in this plot, defin¬ 
ing clearly distinct sequences. As in Figs. |8] and HI Sc 
galaxies constitute a well-defined transition type in this 
diagram, in the sense that earlier-type galaxies are up to 
a factor of ^ 10 more massive than later-type ones for 
the same outermost size i? 25 . 5 - 

Isophotal radii represent a robust metric of global 
galaxy sizes in the nearby universe, but their applicabil¬ 
ity to distant galaxies is hampered by cosmological sur¬ 


face brightness dimming. In order to facilitate the com¬ 
parison of the S^G mass-size relation with observations 
at higher redshits, we have also measured the effective 
radius r^g of our galaxies. In Fig. 1151 we show the trend 
between the effective radius at 3.6 /rm and the total stel¬ 
lar mass of our galaxies, in bins of morphological type. 
As explained in Sect. 15.61 res was measured along the 
semi-major axis of elliptical apertures, so no inclination 
corrections are required. The diagonal lines in this plot 
represent different values of the average stellar mass sur¬ 
face density inside reS, that is, O.bM^,/{irr'^g). Points are 
again color-coded according to their extinction-corrected 
sSFR. 

The main difference with respect to the mass-size re¬ 
lation based on isophotal radii (Fig. [13 is an increased 
spread in M* at fixed size, because a higher concentra¬ 
tion index will significantly shrink r^g, but not so much 
the outer isophotal radius. As a result, early-type galax¬ 
ies can be up to two orders of magnitude more massive 
than late-type ones with the same reff. 

As a comparison, we have overplotted the mass-size 
relations for early - and late-type galaxies found by 
iLange et al.1 ()2014 ) at 0.01 < z < 0.1, based on data 
from the Galaxy and Mass Assembly survey (GAMA, 
[Driver et al.ll201ll) . These authors relied on different cri¬ 
teria to separate early-type galaxies from late-type ones: 
Sersic index, colors and visual identification; here we plot 
the latter. We show their mass-size relation in the r- 
band (where their sample selection was performed), but 
shifting Teff down by 0.076 dex to account for the slightly 
smaller sizes of galaxies at 3.6 p,m. This offset was de- 
rived from the em pirical fits of red vs. wavelength in 
[Lange et al.l ([2014D . 

We find an excellent agreement between the S^G mass- 
size relation and the GAMA one. Massive early-type 
galaxies exhibit a roughly constant average stellar sur¬ 
face density of 10® Mq kpc“^, whereas in late-type galax¬ 
ies the surface density decreases monotonically from 
10® Mq kpc~^ in massive disks to 10^ Mq kpc~^ in low- 
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Fig. 15.— Effective radius along the semi-major axis at 3.6/.tm as a func tion o f the total stellar mass, grouped by morphological type. 
Data-points are color-coded based on the sSFR of each galaxy (see Sect. 16.1.11 for details). The dashed diagonal lines mark constant 
values of the average stell ar mass surface den sity inside rgff. The red and blue contours delineate the mass-size relation of early- and 
late-type galaxies found_b^JLangee^^ <120141) in the GAMA survey. The orange and green curves correspond to the early- and late-type 
distributions found bv IShen e^lTi|2()Q^ ) on SDSS data. The latter were derived from circular apertures; the vertical arrows in each panel 
show how our ellipse-based measurements would shift downwards had we used circular apertures, given the median ellipticity in each bin 
of morphological type. 


mass disks (|Kauffmann et al'1l2003D . The transition from 
the early- to the late-type sequence occurs at a specific 
SFR of ~ 10“^^yr“^. In agreement with our findings 
when discussing the spread in concentration at fixed stel¬ 
lar mass, here we also note that the early-type mass-size 
sequence is not entirely populated by red and old sys¬ 
tems: bars, mergers and XUV-disks can lead to relatively 
high sSFR in galaxies that have otherwise high infrared 
concentration indices and stellar mass surface densities 
(see Fig. [T^. 

A similar mass -size relation was previously found by 
iShen et ^ (|2003D . using SDSS data of nearby galaxies as 
well. In Fig. [12] we have overplotted their trends (±Icr) 
for early- and late-type galaxies, which were classified 
as such based on their concentration and Sersic indices. 
It is worth noting that IShen et HI (|2003n used circular 
apertures instead of elliptical ones, an this typically de¬ 
creases the effective radius as Teffcirc = where 

6/a is the axial ratio. We computed the median offset 


that our data-points would undergo had we used circu¬ 
lar apertures, based on the distribution of axial ratios in 
each bin of morphological type. These offsets are shown 
as vertical arrows in Fig. [15] . Offsets are negligible in 
early-type, spheroid-dominated galaxies, which tend to 
be round regardless of the observing angle. However, in¬ 
clination does play a role in late-type, disk-dominated 
galaxies, where res can decrease by ^ 0.1 — 0.2 dex when 
circular apertures are used. 

6.2.1. Direct size comparisons with previous studies 

Galaxy size measurements, and effective radii in par¬ 
ticular, can be biased by the image depth, the observed 
wavelength (due to radial color gradients), and the fit¬ 
ting methodology. In order to assess the impact of these 
factors in determining accurate galaxy sizes, here we per¬ 
form a direct galaxy-to-galaxy comparison of the S'^G 
effective radii with published measurements from other 
surveys that overlap with ours. 
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Fig. 16.— Comparison between the 2MASS XSC effective radii 
and the ones. For each galaxy, the 2MASS radius is the av¬ 
erage of the effective radii in J, H and Ks^ The S'^G radii are 
measured at 3.6/xm. The color code reflects the average surface 
brightness inside the effective radius at 3.6 pm, without correcting 
for inclination, in AB mag arcsec"^. The inset plot shows the loga- 
ritmic ratio of the 2MASS and S'^G effective radii for galaxies with 
{M)efF3.6 < 21 (blue histogram) and > 21 (red histogram). 

We first cross-matched our sam ple with the 2 MAS S 
Extended Source Catalog (XSC; IJarrett et al.l 120001) . 
which yielded 1687 galaxies in common (unmatched 
galaxies typically have ms.e > 13 — 14 AB mag). The 
XSC provides surface photometry carried out in essen¬ 
tially the same way as in S^C: first, ID radial profiles 
are measured using concentric elliptical apertures; sec¬ 
ond, total magnitudes are obtained by fitting the outer 
part of the growth curve; finally, Teff is derived by inte¬ 
grating on the growth curve until reaching the half-light 
value. The XSC quotes effective radii along the major 
axis in J, H and Ks, but the values differ by less than 
^ 5 — 10% among the three bands, so for each galaxy we 
simply averaged the three values. 

In Fig.[T6]we compare the S^C and XSC effective radii. 
Civen the similar methodologies and wavelengths of both 
surveys one would expect an excellent agreement between 
both sets of measurements. While this is generally the 
case, the XSC sizes are considerably affected by the shal¬ 
lowness of the 2MASS images. To illustrate this, we 
have color-coded the data-points according to (/i)eff3.6i 
the average surface brightness inside the half-light el¬ 
liptical aperture at 3.6/im, uncorrected for inclination. 
For bright galaxies with an effective surface brightness 
{fAeSs.e < 21 AB mag arcsec”^ the 2MASS radii agree 
with the S^C ones with a Itr scatter of ±20%. However, 
galaxies with fainter surface brightness levels appear to 
be several times smaller in 2MASS. This is further high¬ 
lighted in Fig.dTl where we compare the S^C and 2MASS 
images of low surface brightness galaxies whose effective 
radius is severely underestimated in 2MASS. This effect 
must be also taken into account at high redshift, where 
cosmological surface brightness dimming hampers the de¬ 
tection of galaxy outskirts. 

We can also compare our size measurements with 



Fig. 17.— Comparison of S^G 3.6/rm images (top) and 2MASS 
J ones (bottom) for some representative galaxies with (/r)eff 3.6 > 
21 AB mag arcsec”^. Images of the same galaxies are displayed at 
the same spatial scale. The cyan ellipses show the respective half- 
light apertures. 

those obtained from SDSS data using parametric fits, 
in order to understand how many components these fits 
must have in order to accurately recover the global non- 
parametric sizes of galaxies. The SDSS pipeline performs 
simple exponential or de Vaucouleurs fits to the light 
prohles of galaxies, but as we will show below these ba¬ 
sic fits provide only rough (and potentially biased) esti¬ 
mates of galaxy sizes. More sophisticated fits, including 
Sersic profiles and/or multi-component fits, have been 
widely used by other authors on SDSS data (see , e.g. , 
Blanton et al.l boO^ iCadottil l2009t Simard et al.l 120111 : 

Kelvin et al.ll2012l: iLackner &: Cuiml 20121). 

Here we compare our ef f ective radii with those mea¬ 
sured by ILackner fc CunnI (|2012D on a sample of SDSS 
nearby galaxies, which includes 124 S'^C objects. These 
authors fitted the r-band image of each galaxy with five 
different 2D models: a single exponential model, a single 
de Vaucouleurs one, a Sersic model, an exponential bulge 
plus an exponential disk, and a de Vaucouleurs bulge plus 
an exponential disk. 

We compare their effective radii along the major axij^ 
with the S^C ones in Fig. [TS] Panels (a) to (e) show 
the five models de scribed above, and panel (f) shows the 
model deemed by ILackner fc CunnI (120121 ) as the most 
representative for each particular galaxy. The r-band 
radii are systematically larger than the 3.6/im ones by 
^ 0.07 — 0.08 dex. This is quantitatively co nsistent with 
the wa velength dependence of reff found bv iLange et al.l 
(|2014D . due to color gradients. After accounting for this 
offset, in term of scatter the “best model” option in panel 
(f) provides the best agreement with our non-parametric 
sizes, with a la scatter of merely ±0.07dex (±16%). 
The bulge±disk models yield a similar accuracy. The 
Sersic fit does a somewhat poorer job, recovering the 
non-parametric sizes within ±22%. For the single expo¬ 
nential fit, the scatter around the 1:1 line is broader and 
more skewed, as such model is only suitable for almost 
bulgeless systems. Finally, the de Vaucouleurs profile 
completely overestimates the sizes except for some com¬ 
pact objects. 

For the bulge + disk models, ILackner fc Guniil II2012I1 do not 
provide global reff for the whole galaxy, only for each individual 
component. We therefore added up the bulge and disk ID profiles 
along their respective major axes and computed a global reff. 
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Fig. 18.— SDSS r-band effective radii from iLackner fc GunnI 1120121) compared to the S^G ones at 3.6 pm. The Lackner radii result 
from fitting each galaxy with various 2D models: (a) single exponential profile; (b) single de Vaucouleurs profile (the full range is not 
shown to avoid unnecessarily rescaling the axes of the other panels); (c) single Sersic profile; (d) exponential bulge + exponential disk; 
(e) de Vaucouleurs bulge + exponential disk; (f) model that best describes each galaxy. Points are color-coded according to the 3.6 pm 
concentration index C 82 . The small histograms show the distribution of log(refrSDSSr/^eff S4 g)- The mean and Icr values of these 
distributions are quoted in dex units as (p, a) in each panel. 


7. ANCILLARY DATA AND ONLINE ACCESS 

The enhanced data products from our different 
pipelines are available through a dedicated webpage at 
the NASA/IPAC Infrared Science Archiv^ED- For each 
galaxy, users can download the PI images and weight- 
maps, the P2 masks, the P3 radial profiles, the P4 GAL- 
FIT decompositions, and the P5 stellar mass maps. 

IRSA also hosts an interactive catalog that contains 
all the global properties measured in P3 for each galaxy: 
central coordinates, background level and noise, outer 
isophotal sizes and shapes, apparent and absolute asymp¬ 
totic magnitudes, stellar masses, and concentration in¬ 
dices. We have also incorporated into this catalog other 
parameters mined from external sources. In particular, 
the catalog contains redshift-independent distances from 
NED, and a variety of ancillary measurements compiled 
in HyperLeda such as optical sizes, magnitudes and col- 
ors, radial velocities, gas content and internal kinemat- 
icso This allows users to define detailed subsamples 
within S^G according to different selection criteria; the 
resulting list of galaxies can be then fed into the IRSA 
query system to retrieve the images, profiles and other 
data for the selected objects. 

8. SUMMARY 

http://irsa.ipac.caltech.edu/data/SPITZER/S4G/ 

Note that our IRSA catalog is not dynamically linked to NED 
or HyperLeda. Therefore, changes in these external sources will not 
be reflected in the values quoted in the IRSA catalog. 


S^G is the largest and most homogeneous inventory 
of the stellar mass and structure in the nearby universe. 
Our 3.6 p,m and 4.5 /tm data peer through dust at the old 
stellar backbone of more than 2300 nearby galaxies of all 
morphological types in different environments. 

In this paper we have described the methodology and 
algorithms of two of the pipelines in our data flow. 
Pipeline 2 creates masks of foreground stars and back¬ 
ground galaxies. Pipeline 3 measures the background 
level and noise, performs surface photometry on the im¬ 
ages, and derives global quantities such as asymptotic 
magnitudes, stellar masses, isophotal sizes and shapes, 
and concentration indices. 

We find a structural transition between disk- 
dominated to bulge-dominated galaxies at stellar masses 
between 10^° and 10^^ Mq, concurrent with a drop in the 
specific SFR from ^ 10“^°'®yr“^ to ^ 10“^^'®yr“^ with 
increasing concentration. However, we find that not all 
galaxies with high stellar concentration indices are qui¬ 
escent spheroidal systems. Bars and mergers can yield 
central concentrations of stellar mass while at the same 
time enhancing the overall SFR of the galaxy. Extended 
UV disks can also increase the SER in the outer disk of 
galaxies that look otherwise compact in the IRAG bands. 

We have also studied the local stellar-mass size rela¬ 
tion at 3.6 /tm, which is consistent with previous results, 
but has the advantage of being little influenced by gra¬ 
dients in dust extinction or stellar age and metallicity. 
Early-type galaxies present an approximately constant 
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stellar mass surface density of lO®M 0 kpc“^. In late- 
type galaxies the surface density decreases with mass 
from 10® Mq kpc“^ to lO’^MQkpc”®. 

Given the large ancillary value of S'^G, we have made 
our enhanced data products available to the community 
through a dedicated webpage at IRSA. This data release 
includes not only the products described in the present 
paper (masks, radial profiles and derived global quanti¬ 
ties), but also other enhanced products from our various 
pipelines, including science-ready mosaics and weight- 
maps, 2D image decompositions, and stellar mass maps. 
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APPENDIX 

PIPELINE 1 TECHNICAL DETAILS 

Th e S^G Pipeline 1 is an evolved version of the IRAC pipeline used for the SINGS Legacy Program (iKennicutt et al.l 
IM)1 . A major change to the code from previous versions is the removal of a step to improve the registration of 
individual frames using cross-correlation, because the Spitzer calibration now uses 2MASS to improve the astrometry. 
Also, sky subtraction has been moved into Pipeline 3. Below we summarize the steps in the pipeline. 

Saturation removal and rough cosmic ray rejection 

The first step is optional and only occurs when the IRAC observations were taken in High Dynamic Range (HDR) 
mode. In this mode a short exposure is taken before the main exposure at each location. For the vast majority of the 
observations in the sample, the main observation is 30 seconds and the HDR exposure is 1.2 seconds. For new galaxies 
that were only observed during the warm mission we did not use HDR mode since very few have bright nuclei. 

If the Basic Calibrated Data (BCD) images have an HDR counterpart and the flux in the pixel is above 10® MJy sr“^, 
the value from the normal BCD is replaced by the HDR one. In addition, we do a rough cosmic ray removal by 
replacing the value of the pixels in the normal BCD image with the value in the HDR pixel when the difference is 
above 3MJysr“^. 

Optical distortion correction and rotation 

Before we can compare the sky level in overlapping BCD images, we must first correct each BCD for optical distortion 
and rotate the images to have the same orientation. This is required to get pixels to all have the same projeted size 
and orientation on the sky. To do this we use the Drizzle routine and create a new optically corrected BCD image for 
each BCD image in the mosaic. Once we have all the images with constant size on the sky we can find the matching 
regions of overlap. 

Determining the background offset between BCD images 

For all pairs of overlapping BCD images where at least a minimum number of pixels (nominally 20,000) overlap 
between the two images we find the background offset. To determine the background offset we find the difference 
in the 20th percentile pixel flux in both images. Using the 20th percentile has several advantages over a mean or 
median. Our goal is to find the change in the background level rather than any changes when pixels receive flux from 
resolved sources. Sub-pixel changes in the telescope position will change the flux level of pixels observing a point 
source substantially with the IRAC undersampled optics. These same sub-pixel changes do not affect the observations 
of the unresolved background. Also, some detector artifacts such as mux bleed and column pull down can affect large 
areas leading to biases in even the median, while using the 20 th percentile value minimizes the bias. 

For some archive observations the original frames had minimal overlap for the mosaic. In those cases we had to 
reduce the required number of pixels in the overlap region to be as low as 5000. 

Solving for the background offset correction 

Once we have t he ba ckground offset between each pair of overlapping BCD images, we used the method of 
iRegan fc Gruendl (|1995li to find the background offset for each image that minimizes the residuals. This method 
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performs a least square minimization of the residuals using a single offset that is applied to each image. 

Due to their small size, many sample galaxies were not observed in IRAC mapping mode but instead were observed 
in dither mode. In this case two (or even three) non-overlapping regions of the sky are observed and there is a need of 
additional information to obtain a solution. Therefore, when the regions of the sky do not overlap we add in an offset 
between a BCD image in each region of the sky. This will force the sky level to be similar between the non-overlapping 
regions of the sky. 

The set of measured offsets of overlapping images does not provide a zero point for the solution, which requires 
having an external paramenter fixed to set the zero point for the relative offsets. To minimize the effect of applying 
the offsets on the background level, we place the additional constraint that the average offset should have a value of 
0. Once we know the correction that needs to be applied, we create a corrected version of each input BCD image. 

Drizzle the corrected images and update FITS header 

The final mosaics are created using the standard STDAS PYRAF task Drizzle. The resuling images have a plate 
scale of 0.75" per pixel and are oriented to have north to the top. We have to make a correction to the output of 
the drizzling task to recover the correct photometry. Drizzle assumes the pixels have units of flux instead of surface 
brightness. We apply a correction to account for the smaller pixel area in the final mosaic. Drizzle also returns a 
weight map that shows the number of seconds of integration time for each final pixel. This image is also included in 
the S'^G archive and is used by subsequent S^G pipelines. 

In the final step we update the FITS header to include a list of the Astronomical Observation Requests (AORs) that 
contributed to the final mosaic and set all the pixels outside of the field of view to have the value of NAN. 



Surface photometry of galaxies 


19 


REFERENCES 


Aaronson, M., Huchra, J., &; Mould, J. 1979, ApJ, 229, 1 
Abraham, R. G., Tanvir, N. R., Santiago, B. X., et al. 1996a, 
MNRAS, 279, L47 

Abraham, R. G., van den Bergh, S., Glazebrook, K., et al. 1996b, 
ApJS, 107, 1 

Abraham, R. G., van den Bergh, S., & Nair, P. 2003, ApJ, 588, 
218 

Baggett, W. E., Baggett, S. M., &; Anderson, K. S. J. 1998, AJ, 
116, 1626 

Barden, M., Rix, H.-W., Somerville, R. S., et al. 2005, ApJ, 635, 
959 

Bell, E. F., Sz de Jong, R. S. 2000, MNRAS, 312, 497 
—. 2001, ApJ, 550, 212 

Bender, R., Doebereiner, S., Sz Moellenhoff, C. 1988, A&AS, 74, 
385 

Bershady, M. A., Jangren, A., & Conselice, C. J. 2000, AJ, 119, 
2645 

Berlin, E., Arnouts, S. 1996, A&AS, 117, 393 
Blanton, M. R., Schlegel, D. J., Strauss, M. A., et al. 2005, AJ, 
129, 2562 

Bosma, A. 1978, PhD thesis, PhD Thesis, Groningen Univ., 

(1978) 

Brinchmann, J., Chariot, S., White, S. D. M., et al. 2004, 
MNRAS, 351, 1151 

Busko, I. C. 1996, in Astronomical Society of the Pacific 
Conference Series, Vol. 101, Astronomical Data Analysis 
Software and Systems V, ed. G. H. Jacoby & J. Barnes, 139 
Calzetti, D. 2001, PASP, 113, 1449 

Caon, N., Capaccioli, M., Sz D’Onofrio, M. 1993, MNRAS, 265, 
1013 

Carter, D. 1978, MNRAS, 182, 797 

Comeron, S., Elmegreen, B. G., Knapen, J. H., et al. 2011, ApJ, 
741, 28 

Comeron, S., Salo, H., Laurikainen, E., et al. 2013, ArXiv 
e-prints, arXiv:1312.0866 
Conselice, C. J. 2003, ApJS, 147, 1 

Conselice, C. J., Bershady, M. A., Dickinson, M., Sz Papovich, C. 
2003, AJ, 126, 1183 

Courteau, S., Dutton, A. A., van den Bosch, F. C., et al. 2007, 
ApJ, 671, 203 

Dale, D. A., Cohen, S. A., Johnson, L. C., et al. 2009, ApJ, 703, 
517 

de Blok, W. J. G., Walter, F., Brinks, E., et al. 2008, AJ, 136, 
2648 

de Jong, R. S. 1996, A&A, 313, 377 

de Vaucouleurs, G. 1977, in Evolution of Galaxies and Stellar 
Populations, ed. B. M. Tinsley &z R. B. Larson, 43—h 
de Vaucouleurs, G., de Vaucouleurs, A., Corwin, Jr., H. G., et al. 
1991, Third Reference Catalogue of Bright Galaxies. Volume I: 
Explanations and references. Volume II: Data for galaxies 
between 0^ and 12^. Volume III: Data for galaxies between 12^ 
and 24^. 

Driver, S. P., Hill, D. T., Kelvin, L. S., et al. 2011, MNRAS, 413, 
971 

Epchtein, N., de Batz, B., Copet, E., et al. 1994, Ap&SS, 217, 3 
Erwin, P. 2005, MNRAS, 364, 283 
Eskew, M., Zaritsky, D., Sz Meidt, S. 2012, AJ, 143, 139 
Eskridge, P. B., Frogel, J. A., Pogge, R. W., et al. 2000, AJ, 119, 
536 

—. 2002, ApJS, 143, 73 

Fall, S. M., Sz Efstathiou, G. 1980, MNRAS, 193, 189 
Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10 
Fernandez Lorenzo, M., Sulentic, J., Verdes-Montenegro, L., Sz 
Argudo-Fernandez, M. 2013, MNRAS, 434, 325 
Franx, M., van Dokkum, P. G., Schreiber, N. M. F., et al. 2008, 
ApJ, 688, 770 

Fruchter, A. S., Sz Hook, R. N. 2002, PASP, 114, 144 
Gadotti, D. A. 2009, MNRAS, 393, 1531 
Gil de Paz, A., Sz Madore, B. F. 2005, ApJS, 156, 345 
Gil de Paz, A., Madore, B. F., Boissier, S., et al. 2005, ApJ, 627, 
L29 

Gil de Paz, A., Boissier, S., Madore, B. F., et al. 2007, ApJS, 173, 
185 

Hauser, M. G., Arendt, R. G., Kelsall, T., et al. 1998, ApJ, 508, 
25 


Ho, L. C., Li, Z.-Y., Barth, A. J., Seigar, M. S., Sz Peng, C. Y. 
2011, ApJS, 197, 21 

Holwerda, B. W., Pirzkal, N., de Blok, W. J. G., et al. 2011, 
MNRAS, 416, 2401 

Holwerda, B. W., Munoz-Mateos, J., Comeron, S., et al. 2013, 
ArXiv e-prints, arXiv: 1309.1444 

Holwerda, B. W., Munoz-Mateos, J.-C., Comeron, S., et al. 2014, 
ApJ, 781, 12 

Jarrett, T. H., Chester, T., Cutri, R., et al. 2000, AJ, 119, 2498 
Jarrett, T. H., Chester, T., Cutri, R., Schneider, S. E., Sz Huchra, 
J. P. 2003, AJ, 125, 525 
Jedrzejewski, R. I. 1987, MNRAS, 226, 747 

Jogee, S., Barazza, F. D., Rix, H.-W., et al. 2004, ApJ, 615, L105 
Kauffmann, G., Heckman, T. M., White, S. D. M., et al. 2003, 
MNRAS, 341, 54 

Kelvin, L. S., Driver, S. P., Robotham, A. S. G., et al. 2012, 
MNRAS, 421, 1007 

Kennicutt, Jr., R. C. 1998, ARA&A, 36, 189 
Kennicutt, Jr., R. C., Armus, L., Bendo, G., et al. 2003, PASP, 
115, 928 

Kent, S. M. 1985, ApJS, 59, 115 

Kim, T., Gadotti, D. A., Sheth, K., et al. 2014, ApJ, 782, 64 
Knapen, J. H., Shlosman, L, Sz Peletier, R. F. 2000, ApJ, 529, 93 
Kormendy, J., Sz Bender, R. 1996, ApJ, 464, L119 
Kormendy, J., Fisher, D. B., Cornell, M. E., Sz Bender, R. 2009, 
ApJS, 182, 216 

Kormendy, J., Sz Kennicutt, Jr., R. C. 2004, ARA&A, 42, 603 
Lackner, C. N., Sz Gunn, J. E. 2012, MNRAS, 421, 2277 
Lange, R., Driver, S. P., Robotham, A. S. G., et al. 2014, ArXiv 
e-prints, arXiv:1411.6355 

Laurikainen, E., Salo, H., Buta, R., Sz Knapen, J. H. 2011, 
MNRAS, 418, 1452 

Li, A., Sz Draine, B. T. 2001, ApJ, 554, 778 

Lotz, J. M., Primack, J., Sz Madau, P. 2004, AJ, 128, 163 

MacArthur, L. A., Courteau, S., Bell, E., Sz Holtzman, J. A. 

2004, ApJS, 152, 175 

MacArthur, L. A., Courteau, S., Sz Holtzman, J. A. 2003, ApJ, 
582, 689 

Marinova, I., Sz Jogee, S. 2007, ApJ, 659, 1176 
Martm-Navarro, I., Bakos, J., Trujillo, L, et al. 2012, MNRAS, 
427, 1102 

Meidt, S. E., Schinnerer, E., Knapen, J. H., et al. 2012, ApJ, 744, 
17 

Meidt, S. E., Schinnerer, E., van de Ven, G., et al. 2014, ApJ, 

788, 144 

Menendez-Delmestre, K., Sheth, K., Schinnerer, E., Jarrett, 

T. H., Sz Scoville, N. Z. 2007, ApJ, 657, 790 
Minchev, I., Famaey, B., Combes, F., et al. 2011, A&A, 527, A147 
Mo, H. J., Mao, S., Sz White, S. D. M. 1998, MNRAS, 295, 319 
Morgan, W. W. 1958, PASP, 70, 364 

Munoz-Mateos, J. C., Boissier, S., Gil de Paz, A., et al. 2011, 

ApJ, 731, 10 

Munoz-Mateos, J. C., Gil de Paz, A., Boissier, S., et al. 2007, 

ApJ, 658, 1006 

Munoz-Mateos, J. C., Gil de Paz, A., Zamorano, J., et al. 2009a, 
ApJ, 703, 1569 

Munoz-Mateos, J. C., Gil de Paz, A., Boissier, S., et al. 2009b, 
ApJ, 701, 1965 

Munoz-Mateos, J. C., Sheth, K., Gil de Paz, A., et al. 2013, ApJ, 
771, 59 

Paturel, G., Petit, C., Prugniel, P., et al. 2003, A&A, 412, 45 
Peletier, R. F., Davies, R. L., Illingworth, G. D., Davis, L. E., Sz 
Cawson, M. 1990, AJ, 100, 1091 
Peng, C. Y., Ho, L. C., Impey, C. D., Sz Rix, H.-W. 2002, AJ, 

124, 266 

—. 2010, AJ, 139, 2097 

Querejeta, M., Meidt, S. E., Schinnerer, E., et al. 2014, ArXiv 
e-prints, arXiv:1410.0009 

Regan, M. W., Sz Gruendl, R. A. 1995, in Astronomical Society of 
the Pacific Conference Series, Vol. 77, Astronomical Data 
Analysis Software and Systems IV, ed. R. A. Shaw, H. E. 

Payne, &; J. J. E. Hayes, 335 

Regan, M. W., Sheth, K., Sz Vogel, S. N. 1999, ApJ, 526, 97 
Roskar, R., Debattista, V. P., Stinson, G. S., et al. 2008, ApJ, 

675, L65 


20 


Munoz-Mateos et al. 


Rubin, V. C., Thonnard, N., & Ford, Jr., W. K. 1978, ApJ, 225, 
L107 

Ryden, B. 1992, ApJ, 396, 445 
Ryden, B. S. 2004, ApJ, 601, 214 
—. 2006, ApJ, 641, 773 

Salo, H., Laurikainen, E., Laine, J., et al. 2015, ArXiv e-prints, 
arXiv: 1503.06550 

Sanchez-Blazquez, P., Courty, S., Gibson, B. K., & Brook, C. B. 
2009, MNRAS, 398, 591 

Sandage, A., Freeman, K. C., & Stokes, N. R. 1970, ApJ, 160, 831 
Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 
525 

Scodeggio, M., Gavazzi, G., Franzetti, P., et al. 2002, A&A, 384, 
812 

Shen, S., Mo, H. J., White, S. D. M., et al. 2003, MNRAS, 343, 
978 

Sheth, K., Regan, M. W., Vogel, S. N., & Teuben, P. J. 2000, 
ApJ, 532, 221 

Sheth, K., Vogel, S. N., Regan, M. W., Thornley, M. D., & 
Teuben, P. J. 2005, ApJ, 632, 217 
Sheth, K., Elmegreen, D. M., Elmegreen, B. G., et al. 2008, ApJ, 
675, 1141 

Sheth, K., Regan, M., Hinz, J. L., et al. 2010, PASP, 122, 1397 
Simard, L., Mendel, J. T., Patton, D. R., Ellison, S. L., & 
McConnachie, A. W. 2011, ApJS, 196, 11 
Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 
1163 

Sofue, Y., & Rubin, V. 2001, ARA&A, 39, 137 
Sorce, J. G., Courtois, H. M., & Tully, R. B. 2012, AJ, 144, 133 
Taylor, V. A., Jansen, R. A., Windhorst, R. A., Odewahn, S. C., 
& Hibbard, J. E. 2005, ApJ, 630, 784 


Taylor-Mager, V. A., Conselice, C. J., Windhorst, R. A., &: 
Jansen, R. A. 2007, ApJ, 659, 162 

Thilker, D. A., Bianchi, L., Boissier, S., et al. 2005, ApJ, 619, L79 

Thilker, D. A., Boissier, S., Bianchi, L., et al. 2007, ApJS, 173, 
572 

Trujillo, I., Rudnick, G., Rix, H.-W., et al. 2004, ApJ, 604, 521 

Trujillo, I., Forster Schreiber, N. M., Rudnick, G., et al. 2006, 
ApJ, 650, 18 

Tully, R. B., & Fisher, J. R. 1977, A&A, 54, 661 

van der Wei, A., Franx, M., van Dokkum, P. G., et al. 2014, ApJ, 
788, 28 

Verheijen, M. A. W. 2001, ApJ, 563, 694 

Wang, J., Kauffmann, G., Overzier, R., et al. 2011, MNRAS, 412, 
1081 

Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, ApJS, 154, 
1 

Whyte, L. F., Abraham, R. G., Merrifield, M. R., et al. 2002, 
MNRAS, 336, 1281 

Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, 
AJ, 140, 1868 

York, D. G., Adelman, J., Anderson, Jr., J. E., et al. 2000, AJ, 
120, 1579 

Younger, J. D., Cox, T. J., Seth, A. C., & Hernquist, L. 2007, 
ApJ, 670, 269 

Zaritsky, D., Salo, H., Laurikainen, E., et al. 2013, ApJ, 772, 135 

Zaritsky, D., Courtois, H., Munoz-Mateos, J.-C., et al. 2014, AJ, 
147, 134 


